Band structure renormalization and weak pseudogap behavior in Na .33CoO2: 
Fluctuation exchange study based on a single band model 
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Based on a single band Hubbard model and the fluctuation exchange approximation, the effective 
mass and the energy band renormalization in Nao.33Co02 is elaborated. The renormalization is 
observed to exhibit certain kind of anisotropy, which agrees qualitatively with the angle-resolved 
photoemission spectroscopy (ARPES) measurements. Moreover, the spectral function and density 
of states (DOS) in the normal state are calculated, with a weak pseudogap behavior being seen, 
which is explained as a result of the strong Coulomb correlations. Our results suggest that the large 
Fermi surface (FS) associated with the ai g band plays likely a central role in the charge dynamics. 
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The layered oxide NazCoC^ has attracted much at- 
tention due to its possible connection to the high-T c su- 
perconductivity since the discovery of superconductivity 
in itfii^ This material consists of two-dimensional C0O2 
layers, where Co's form a triangular lattice, making the 
Na x Co02 a possible realization of Anderson's resonating 
valence bond (RVB) stated By hydration, it becomes a 
superconductor with T c w 5K in Nao.337CoO2d.3H2O 1 . 
The unhydrated compound Na 2: Co02 exhibits a rich 
phase diagram: a paramagnetic metal (x < 0.5), a 
charge-ordered insulator (x = 0.5), a "Curie- Weiss 
metal" (x ~ 0.7), and a magnetic order state (x ^ 0.75). 4 

Experimentally, it is indicated that Na K Co02 seems to 
be strongly electronic correlated^^ Angle-resolved pho- 
toemission spectroscopy (ARPES) measurements show a 
strong mass renormalization (The effective mass is about 
5-10 times larger than the bare mass.) for x = 0.6Z De- 
creasing the Na concentration to x = 0.3, the cobaltates 
appear to be less electron-correlated with a weaker but 
still apparent energy band renormalization factor (~ 2). 8 
For the hydrated compounds, the effective mass is esti- 
mated to be ~ 2 — 3^ In addition, recent experiments 
reported that Nao.33CoO2-l.3H2O display certain pseu- 
dogap behaviors such as the decreasing of the Knight 
shift 10 and the density of states (DOS) at the Fermi level 
below 300K4 1 - Optical spectroscopy measurements for 
Na .2CoO 2 and Na .5CoO 2 also suggest the incipient for- 
mation of pseudogap. 12 However, unlike the pseudogap 
effect observed in the high-T c cuprates, the pseudogap be- 
havior is rather weak. The renormalization of the energy 
band and the pseudogap formation are directly related to 
the electronic structure, such as the quasiparticle spec- 
tral function, the quasiparticle dispersion, and the Fermi 
surface (FS) topology. Therefore it is of importance and 
significance to study the normal state quasiparticle dy- 
namics. 

It is believed that the topology of the FS plays an 
important role in the unconventional superconductivity. 
The local density approximation (LDA) calculation s 13 ' 14 



predict a large FS associated with the a\ g band and 
six pockets associated with the e' g band. However, in- 
triguingly, the small pockets have not been observed in 
the ARPES measurements*^^ The unexpected incon- 
sistency of the topology of the FS between the LDA re- 
sults and the ARPES measurements has aroused much 
controversy. One possibility might be due to the surface 
effect in the ARPES measurements. From another view- 
point, Zhou et al. suggested that the strong Coulomb in- 
teractions may induce significant renormalization of the 
band structure which pulls the e' g band down from the 
Fermi energyii 6 Furthermore, the e' g band is suggested 
to be relevant to the spin-triplet superconductivity in 
some theoretical and numerical studies! 17 ' 18 On the other 
hand, the possibility of the spin-singlet superconductivity 
is also suggested according to recent Knight-shift mea- 
surements ^2iS£i Taking the above considerations into ac- 
count, a single band model focusing on the a\ g band may 
be a starting point as a minimal model of Na x Co02- 

In this paper, we study the normal state electronic 
structure in Nao.33Co02 on the basis of the single band 
Hubbard model and the fluctuation exchange (FLEX) ap- 
proximation. Several experimental features, such as the 
band renormalization and the weak pseudogap behavior, 
are well reproduced, which suggests an important role of 
the holelike Fermi surface centered around the T point 
in the quasiparticle dynamics. We start with the two- 
dimensional single band Hubbard model given by, 



H^2n la (1) 



where tjj denotes the hopping term (in the following we 
will use t\,t2, ts and £4 to denote the hoppings along the 
nearest- to the fourth nearest-neighbors), U the on-site 
Coulomb repulsion and fi the chemical potential. To re- 
produce the large FS around the T point, we set the 
parameters (t±, <2, *3, £4) of the bare dispersion ek to 
be (-1, 0, 0, 0.2), where t\ is set to be the unit hereafter. 
According to the numerical result for the bandwidth ; 13 ' 21 




FIG. 1: The Fermi surface of Na .33CoO 2 for (ti,t 2 ,t 3 U)={- 
1,0,0,0.2). The grey area denotes the unoccupied electron 
states. Point A near the Fermi surface indicates the k-point 
where m* is calculated. 



one can get t\ = 130 meV. We note that the main feature 
of the bare dispersion is similar to that of the band with 
an aig character in the LDA calculation by Lee et al.^ 
The FLEX approximation is employed in our study, 
which has been applied to the two dimensional Hubbard 
model in literatures i 22 i 23 ' 24 i 25 As a self-consistent approx- 
imation, the FLEX approximation solves the Dyson's 
equation with a primarily RPA+ladder type effective in- 
teraction self-consistently. Based on the scenario of the 
FLEX approximation, the self-energy is given by, 



E(fc) 
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where 
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The spin susceptibility is 

Xs(q) = 



x(q) 



1 - Ux(q) ' 
and the charge susceptibility is 
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The electron Green's function is give by, 

G(k) = [iUn-Ck-Sik)]- 1 . 
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In the above equations, k = (k,iu> n ) and q = (q, iv n ) 
are used, T is the temperature. These equations are 



FIG. 2: The ratio of the effective mass to the bare mass m*/m 
at point A as indicated in Fig. 1 versus the on-site Coulomb 
repulsion U at T — 0.1. 



solved self-consistently, where TV = 64 x 64 k-point 
meshes and up to 2048 Matsubara frequencies u> n — 
(2n + 1)ttT are taken. The electron density is deter- 
mined by the chemical potential fi from the equation 

<n>=l-fE t (?(t). 

As an important physical quantity related closely to 
the strength of the electronic correlations, the effective 
mass m* is given by 



m „. m d „ 
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is the renormalization constant. The self-energy with real 
frequencies E(k, w) is obtained with the Pade approxi- 
mants, 2 - which analytically continue S(k, iui) from the 
Matsubara frequencies to the real-frequency axis. The 
ratio of the effective mass to the bare mass m*/m with 
different on-site Coulomb repulsions U near the FS (point 
A in Fig. 1) is plotted in Fig. 2, where m* jm increases 
monotonously with the increasing of U. It is noticeable 
that when U is less than 4, the enhancement of the effec- 
tive mass is inapparent (< 1.2). With a strong Coulomb 
repulsion (U > 6 w 0.8 eV), the effective mass is en- 
hanced significantly. Comparing the calculated effective 
mass with the ARPES measurements (~ 2 for x ~ 0.3), 
we set U ~ 6 in the present simple model to reflect 
the strong correlation effect in Nao.33Co02^ The strong 
electronic correlations are evident from Fig. 3 (a), where 
the sharp quasiparticle peak is suppressed with strong 
Coulomb repulsion. For U — 2 and T — 0.1, the spec- 
tral function shows a sharp peak near the Fermi energy. 
Increasing the Coulomb repulsion to U = 6, the spec- 
tral function is broadened and the spectral weight at the 
Fermi level is reduced greatly. 

For U = 6 we now investigate the renormalization of 
the quasiparticle energy band. To do it, we need to eval- 
uate the spectral function of quasiparticles defined as 
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FIG. 3: (a) The effect of the on-site Coulomb repulsion on 
the spectral function at point A of Fig. 1. Solid line: (7 = 6, 
T = 0.1. Dashed line: U = 2, T = 0.1. (b) Solid line: the 
renormalized band dispersion for U — 6 and T — 0.1. Dashed 
line: the bare band dispersion. Dotted line: the Fermi energy. 



where G(k, w) is the dressed Green's function in real fre- 
quecy, which is analytically continued from the dressed 
Green's function in imaginary-frequecy via the Pade ap- 
proximants. Then, the renormalized band dispersion is 
determined by the position of the peak of A(k, u>). The 
calculated result at T = 0.1 is shown in Fig. 3 (b), 
where the solid and the dotted lines denote the renor- 
malized and the bare energy bands, respectively. One 
can see that the bandwidth is compressed obviously due 
to the strong Coulomb interactions. From Fig. 3 (b) we 
get the bandwidth of the bare band dispersion (7.7) and 
the renormalized bandwidth (5.0), which gives a band 
renormalization factor 1.54. This is consistent with the 
ARPES measurements for x = 0.3. 8 Moreover, we find 
that the band renormalization shows an anisotropy along 
the T-K and T-M direction, namely the renormalization 
along the T-K direction is stronger than that along the 
T-M direction. The similar anisotropy has been found 
in the ARPES experiments for x = 0.6f£ though the ex- 
perimental result seems stronger than what we disclosed 
here. We note that the compound Na.;;Co02 with x = 0.6 
is in fact in the range of the Curie- Weiss metal, in which 
the correlation is stronger than that in the paramagnetic 
metal with x — 0.3. So, a weaker anisotropy in the renor- 
malization is expected. This anisotropy originates from 
the nesting of the Fermi surface along the T-K direction 
as shown in Fig. 1. 

We now turn to address the weak pseudogap behavior. 
This is manifested in the suppression of the density of 
states at the Fermi level. We present the u> dependence 
of the density of states at different temperatures in Fig. 
4. With the decrease of temperature from T = 0.7 (dash 
dotted) through T = 0.5 (dotted) and T = 0.3 (dashed) 
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FIG. 4: The total density of states (DOS) versus energy for 
(7 = 6. Solid line: T = 0.1, dashed line: T = 0.3, dotted line: 
T = 0.5, and dash-dotted line: T = 0.7. Inset: DOS plots for 
U = 2 at T = 0.1, 0.3, 0.5, 0.7. 



to T = 0.1 (solid), a weak suppression of the density 
of states near the Fermi energy is evident, suggesting 
the opening of a weak pseudogap. The weak pseudo- 
gap behavior is also manifested in the quasiparticle spec- 
tral function A(k, ui) which measures the probability to 
find the quasiparticle at momentum k and frequency to. 
As shown in Fig. 5 (a), with the decrease of tempera- 
ture (from the dashed line to the solid line), the spectral 
weight is transferred away from the region around lu = 0, 
producing a weak secondary maximum at uj < Ef- Thus, 
a weak pseudogap is formed. 

This weak pseudogap behavior is a consequence of 
strong Coulomb repulsion. To show this, we also present 
the results with a smaller U = 2 at T = 0.1,0.3,0.5,0.7 
in the insert of Fig. 4. Different from the U = 6 case, the 
density of states is not suppressed with the decreasing of 
temperature at all(in fact there is no appreciable differ- 
ence between them, so only one line can be seen from the 
figure). For a more detailed discussion, we also present 
the spectral function for U = 2 (dotted), U — 4 (dash 
dotted) and U = 6 (solid) in Fig. 5 (a), where the sec- 
ondary maximum shows up gradually with the increase of 
the Coulomb correlations. Note that since the secondary 
maximum disappears for U — 2, no pseudogap is ex- 
pected for this case and what for U < 2. We refer to the 
depression of the spectral weight to be a weak pseudogap 
here is based on the observation that the spectral func- 
tion consists of a peak and a weak secondary maximum 
(not a peak). In fact, the real part of the denominator 
of the Green's function has only one pole, which can be 
seen from the real part of the self-energy shown in Fig. 5 
(b). It is the extremum of the imaginary part of the self 
energy around oj = — 4 that produces the weak maximum 
of the spectral function, as shown in Fig. 5 (c). In the 
case of a smaller Coulomb repulsion U, such as U = 2, 
a usual Fermi-liquid form of the self-energy is preserved 
(The dotted lines in Fig. 5 (b) and (c)), so it gives a 
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FIG. 5: (a) The energy dependence of the spectral function 
at the momentum indicated as point B in Fig. 1. Below the 
Fermi energy, there is a small peak in the spectral function 
for U = 6 and T = 0.1(solid line). Dashed line: (7 = 6 
and T = 0.7, dash-dotted line: U = 4, T = 0.1, dotted 
line: U = 2, T = 0.1. The arrow points to the location 
of the secondary maximum. Inset: Zoom-in view of the solid 
rectangle, (b) and (c): The real and the imaginary part of the 
self-energy at T = 0.1 for U = 6 (solid) and U = 2 (dotted). 
The arrow indicates the position of the extrema in Im E(k, u>). 



well defined quasiparticle peak in the spectral function 
(The dotted line in Fig. 5 (a)). This secondary max- 
imum might suggest a "shadow band",— which occurs 
when a short-range spin correlation is developed with 
the increase of Coulomb repulsion. Therefore, the weak 
pscudogap behavior reported here may be due to the spin 
fluctuation in the strong correlated regime.— 

A similar weak pseudogap behavior was reported by 
Yada et al£Q based on a multi-orbital model with the 
absence of small hole pockets. Our results based on the 
single band model provide further support for the sinking 
down of the small hole pockets. In our opinion, the agree- 
ment between the two models suggests that the multi- 
orbital effect may play a minor role in the mechanism of 
the pseudogap formation. 

In summary, we have studied the quasiparticle 
band renormalization and the pesudogap behavior in 
Nao.33Co02 based on the single band Hubbard model 
in a two-dimensional triangle lattice. The renormaliza- 
tion of the band structure and its anisotropy of the T-K 
and T-M direction have been elaborated. The estimated 
effective mass is consistent with the ARPES measure- 
ments. In the meantime, a weak pseudogap behavior is 
reproduced, which is explained as the result of the strong 
spin fluctuations. Our results are qualitatively consistent 
with experiments as well as the theoretical calculations 
based on a multi-orbital model. 
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